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Abstract 

A detailed study of the solar neutrino vacuum oscillation is made taking into account three 
neutrino flavours and seasonal effect. A set of Sm^ x sin 2 (2cj) regions is calculated for a range 
of the parameter sin 2 (20) from to 1, with and without the inclusion of the recoil-electron 
spectrum in the rates. The averaged survival probabilities for v e as a function of the energy 
are obtained, what reveal that solutions with values of sin 2 (20) > 0.50 and the maximum 
exclusion of 7 Be neutrinos and minimum exclusion of pp neutrinos give a better explanation for 
the suppression rates of all detectors. 
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1 Introduction 



For three decades the deficit on the observed flux of solar neutrinos compared to that one expected 
from Standard Solar Models (SSM) has puzzled the scientists and originated several explanations 
for this problem. These solutions are divided into astrophysical and beyond the standard model 0. 
Nowadays, the astrophysical solutions have been highly disfavoured by the helioseismological data 
which are well explained by standard solar models, however, solutions with new physics, in special 
the neutrino oscillations, that can be in vacuum or in matter (by MSW effect) ||, gives a good 
fit to the solar neutrino data. 

In the present work we will focus attention only on vacuum oscillations, which occur in the 
neutrino way from the Sun to Earth. The evidence of the existence of vacuum oscillations come from 
the observation of the seasonal effect. Because the Earth orbit is not circular, the variations in the 
Earth-Sun distance produce a difference in the oscillation probability what can be very important for 
neutrino line sources, in special, the 7 Be (0.861 MeV) and is to be confirmed by the Borexino detector 
|4[] . In our study we will deal with two and three neutrino flavours to try to explain the solar neutrino 
problem. We will vary the parameter sin 2 (20) that is concerned with the z/ e ^z/ T oscillations in order 
to discover which scenario (two or three generations) gives a better description for the experimental 
data. This analysis will be made by the calculation of the 8rn^ 1 x sin 2 (2a;) regions considering first 
only the rates of the detectors without the spectral distortion, caused by the recoil-electron spectrum 
in Kamiokande and Super-Kamiokande, and after including these spectrum within the calculation. 
An investigation in the v e survival probability will be made, what will reveal the behavior of the 
neutrino function of the energy. 

In section 2 we give the relations needed to calculate the mass regions and mixing angles allowed 
by the experiments. In section 3 we present the v e survival probability as well the and v T conversion 
probabilities. In section 4 we analise the obtained data to describe the probabilities as a function of 
the energy for a set of values of sin 2 (20). In section 5 we present the conclusions of this work. 



2 Suppression Rates 

The first procedure is to calculate the suppression rates of the experiments using neutrino os- 
cillations and compare with the observed ones to search for the values of the parameters Srn^i and 
sin 2 (2a;) that fit the five experiments (Homestake ||, GALLEX ||, SAGE 0, Kamiokande || and 
Super-Kamiokande 0) together. The suppression rate, to the experimental case, is obtained by di- 
viding the capture rate of the detector to that one expected from the SSM and theoretically, dividing 
the capture rate calculated using neutrino oscillations to the SSM one. We will make this calculation 



using the data from the BP98 Standard Solar Model [10| with 99% C.L. 

The Solar Neutrino Problem is evidenced when we compare the theoretical capture rates with 
those obtained experimentally. When this comparison is made we discover that these rates never 
match, giving values varying from 33% for Homestake to 60% for GALLEX. See table 1. 

The solar neutrino suppression rate using neutrino oscillation is given by 
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Table 1: Solar neutrino data and BP98 SSM predictions [10]. The second and third columns show 
the measured and calculated capture rates for each detector. The values are given in SNU except 
for Kamiokande and Super-Kamiokande, where the parameter is the 8 B neutrino flux, in units of 
10 6 cm~ 2 s _1 . The suppression rates are set on the fourth column for 1 a errors. 



Experiment 


Result 


BP98 SSM 


Result/SSM 


Homestake []5J 
GALLEX @ 
SAGE 
Kamiokande || 
Super-Kamiokande || 


2.56 ±0.16 ±0.16 
77.5 ± 6.2+^7 
66.6t£? 

2.80 ±0.19 ±0.33 
2.44 ± 0.05j£o7 


7 7+1.2 
' • ' -1.0 

129l|° 
129l|° 
5-lSlo.? 

5 .i5ir 7 


0.33 ±0.03 
0.60 ±0.06 
0.52 ±0.06 
0.54 ±0.07 
0.47 ±0.02 



-'SSM 

where 



R=^ (1) 



max 



Sosc = ^2®i r]i(v v )(T{v v )P VeVR {v v )dE v (2) 

2 ^ E miii 

and S SSM is the same expression with P VaV<s (m v ) = 1 

The equation (§) is available for Homestake, GALLEX and SAGE experiments, where % runs for 
each neutrino source, $j is the total neutrino flux from the source i, ^(e v ) is its normalized neutrino 
energy spectrum, <T{e v ) is the cross section for the experiment considered and P UeUe (E u ) is the electron 
neutrino survival probability. The neutrino sources that we consider for these experiments are pp, 
7 Be, S B, pep, 13 N and 15 0, except for Homestake that is not sensitive to pp neutrinos. 

For Super-Kamiokande we consider the 8 B neutrino source only and we have to take into account 
the fact that this experiment is also sensitive to the v^, e and v T , e scattering, so S osc becomes 

v e ( E y)C"i/„ ( E y) 

dE v (3) 

-Win 

since o v ^ v) = o Vt & v ) and P u ^ e ^ v ) + P Vt v b <p v ) = 1 - Pv e v e (Pu) the equation (|) is simplified to 

r\(v v ) \ P VeVe <v v )O u ^ v ) + (1 - P Ve . Ve .{* v ))<J v „& v ) \ dE v (4) 

-'min 

For the Kamiokande experiment only the 8 B source is considered again but the sensitivity to v T 
is neglected. 

As Super-Kamiokande and Kamiokande are neutrino-electron scattering detectors, different from 
the other three which use neutrino absorption, we can also take into account the energy resolution 
in order to observe the influence of the spectral distortion caused by the recoil-electrons in the 
suppression rates, giving the following expressions for S osc and S SSM 
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rT max m c ^ fT'xnax I _ d<J I/e (T',E u ) 

,(E V )- 



S osc = iV $ / dT^ ( E„) / dT / dT' R{i' ,t) P UeUi 



/"Emax /"Tmax m e /"T'max fxfT (T' F "i 

S SSM = iV $ / dT^ (El/) / dT / dT'Rw v t- v) (6) 

JE min JT min -m c JO al 

The same procedure used to simplify the eq.(|3]) is taken here, what leaves eq.(|) as 

/•Emax fTmax-me /'T'max | (j{J frnj p ^ 

S osc = iV $ / dE uV(Ev) / dT / dT / J R ( T',T) P^ce,)- ' ' 

** E m ; n </ T m ; n — m c J I 



dT 

The energy resolution is described by 

1 [ (T'-T + 5) 2 

R(t',t) = -=exp 



dT 

, p d^ (T', Et/ ) | 



where T and T are the measured and true electron kinetic energy, respectively, and the limit T^ 
in the third integral of the eq.(^) is the maximum kinetic energy that an electron can achieve given 
the neutrino energy E u , T max = E v j(l + m e /2E u ). For more details see [JTTJ, 



The limits of energy used to calculate these rates are not the same the neutrinos have in the Sun, 
so we have to calculate the normalized detector sensitivity to know these limits, as follows 

J e 0(v v )r}i(E v )ah, v 

We can notice in the figure 1 the existence of regions where the spectral sensitivity is zero. 

3 Electron Neutrino Survival Probability 

For three neutrino flavours the v e survival probability is written as W2 



JW,> = 1 - ^Xsin 2 (^ ) - ClSl^i 5 ^) - ^sin 2 (^) (10) 

where the terms C^, Sf w , S" 2 , and Sf^ are short forms for cos 4 (0), sin 2 (2u;), sin 2 (c<j), cos 2 (a;) and 
sin 2 (20), respectively. 

In this work we will consider the following mass hierarchy [TJ] 



ml << m2 << m3 (11) 
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what leads us to Sm^ » 5m 21 and 5m\ 2 ~ Sm^, since 5m| 2 = Sm^ — 6m 21 , thus we can simplify 
the equation (|H]) to the form 



W.) = 1 - CjfiL sin 2 (^) - S 2 2 , sin 2 (^l) (12) 



4hcE ) \sin( 4M 

Now we have two choices 

i) if 5m 21 < 10 _11 ey 2 then the term sin 2 ( ) tends to be very small and we fall in the case of two 
generations with <fi as the mixing angle and Sm^ as the mass parameter; 

ii) if 5m 21 ~ lO -10 ^ 2 , for the range of energy considered the term sin 2 ( S ^ 31 ^ ) can be averaged to 
1/2, giving the final equation for the v e survival probability, which applies for two and three neutrino 
oscillations 

Pu M = 1 - \SI, - Cpl sm 2 (^f ) (13) 

The introduction of the seasonal effect causes a modification in the probability since the distance 
x varies with the time on the following way [14| 

X(t) = x (l - e cos(-^)) (14) 

where x is the mean Earth-Sun distance (~ 1.49 . 10 n m), e is the elipticity of the orbit and T = 365 
days. Once we are taking a temporal average over the probability the eq.(|13"D turns. 

p 1 X Q2 ^4 C 2 1 A*' • 2 / j^lgo(l-eCOS(g)) \ 

P»»M = 1 " " C * S *>Jjt^t3 L Sm 1 ) dt (15) 

where we integrate in a time interval [ti,^]; which in this work we consider as one year. 
For P„^ e (t,E v ) and P VrUe {^ v ) we have 

^2 C 2 1 [ t2 ■ 2 j Sm 2 21 x (l-e cos( 2 ^)) | 

c ^Tb=t5k sm 1 47^ — — r (16) 



and 

p 1 Q2 1 2 c2 1 T 2 • 2 \ 5m 2 21 x {\ -e cos( 2 ^)) l 
P — ™ = 2 S ** - 4^ (^) L Sm I 4teE ' * (17) 

In figure 2 it is shown the expected maximal and minimal variation of these probabilities as a 
function of sin 2 (20). 
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Table 2: Best fit points for the Sm^x x sin 2 (2o;) plots shown on figure 3. Column 1 is for the seasonal 
effect with the rates only and column 2 is with the inclusion of the recoil-electron spectrum. 



sin 2 (20) 


Aram 


sin 2 (2c;) 


5rri2i 


[eV 2 


) 


Column 1 





9.29 


0.76 


6.22 


10" 


-11 


0.25 


9.95 


0.85 


3.17 


10" 


-10 


0.50 


10.3 


0.79 


3.25 


10" 


-10 


0.75 


10.9 


0.91 


3.50 


10" 


-10 


1 


12.9 


0.81 


3.70 


10" 


-10 


Column 2 





4.42 


0.79 


2.70 


10" 


-10 


0.25 


4.78 


0.72 


2.70 


10" 


-10 


0.50 


5.17 


0.96 


7.88 


10" 


-10 


0.75 


5.49 


0.79 


8.15 


io- 


-10 


1 


6.49 


0.63 


8.15 


io- 


-10 



4 Analysis of the Results 

The figure 3 illustrates the plots of Srn^i x sin 2 (2u;) considering the seasonal effect. We can 
notice an increase in the regions that explain the five experiments together as sin 2 (20) runs from 
to 1. The first column shows the regions without taking into account the recoil-electron spectrum 
of Kamiokande and Super-Kamiokande and the second column gives the regions with the inclusion 
of that spectrum, what causes a decrease in the mass regions in the order of 33% for sin 2 (20) = 



to 23% for sin (20) = 1. The x analysis was made following the references [15] and is shown in 
table 2 and we can observe that the inclusion of the recoil-electron spectrum gives a better fit. 

With the Sm 2 ,! and sin 2 (2a;) points from figure 3 we can obtain the averaged electron neutrino 
survival probabilities as a function of the energy, what are ploted on figure 4. In a first sight we 
observe a considerable decrease of the probability on the energy region dominated by the 8 B neutrino 
and a slight decrease in the 7 Be(0.861MeV) line source, because it lies in a region where there are 
CNO neutrinos also. As the parameter sin 2 (20) grows, the probability range decrease with the 
energy. For sin 2 (20) = 0, 0.30 < P„ e „ e < 0.80 and for sin 2 (20) = 1, 0.39 < P VeVe < 0.42. 



To make a better analysis it is useful to adopt the procedure given by Barger et al, Hl3fl , where 
the suppression rates of the detectors are given by means of the averaged probabilities. So we have 
for the detectors the following expressions, 
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Table 3: Suppression rates for the experiments using the averaged solar neutrino survival probabili- 
ties, with 99% C.L. 



sin 2 (20) 


Homestake 


Gallium 


Kamiokande 


Super-Kam 


rates only 





0.52 ±0.19 


0.56 ±0.08 


0.49 ±0.10 


0.50 ±0.11 


0.25 


0.45 ±0.11 


0.52 ±0.04 


0.46 ±0.11 


0.49 ±0.09 


0.50 


0.42 ±0.07 


0.47 ±0.03 


0.45 ±0.09 


0.48 ±0.07 


0.75 


0.41 ±0.04 


0.43 ±0.02 


0.44 ±0.05 


0.48 ±0.04 


1 


0.41 ±0.01 


0.41 ±0.01 


0.43 ±0.02 


0.48 ±0.01 


rates + recoil-electron spectrum 





0.50 ±0.17 


0.55 ±0.08 


0.52 ±0.13 


0.53 ±0.14 


0.25 


0.45 ±0.12 


0.49 ±0.06 


0.51 ±0.15 


0.52 ±0.12 


0.50 


0.43 ±0.09 


0.45 ±0.03 


0.47 ±0.13 


0.50 ±0.10 


0.75 


0.43 ±0.06 


0.42 ±0.02 


0.44 ±0.08 


0.50 ±0.06 


1 


0.43 ±0.03 


0.42 ±0.01 


0.44 ±0.03 


0.50 ±0.02 



R Hom = 0.76(P UeVe )s B + 0.15{P VeVe )7 Be + 0.09( Pu e u e )cNO,pep 

R Gal = 0.l0(Pu e u e )sB + 0.27(P UeUe )7 Be + 0.09(P UeUe )cNO,pep + 0M(P UeUe ) pp ^ 

P-Karri = (Pu e u e ) « B + 7 (Pu^u e ) &B 
Rs-K = 7 ± ^(Pu e Ue) S B 

where for simplification on the calculus we take into account the fact that a Ve , e (^v) ~ 1<j v e ( E v) m 
Kamiokande and Super-Kamiokande. In table 3 there are the values (99% C.L.) of the suppression 
rates using the probabilities from figure 4 considering the following range of energies: E v < 0A20MeV 
for pp and E v = 0.861MeV^ for 7 Be. For 8 B we used E v > 5.0MeV for the radiochemical detectors 
and E v > 6.5 and 7.5MeV for Super-Kamiokande and Kamiokande, respectively. For CNO neutrinos 
the range is 0.814 < E v < 1.732MeV for Homestake and 0.233 < E v < 1.732MeV for the gallium 
detectors. 

A comparison between table 3 and the experimental data from table 1 shows us that the sup- 
pression rates for Homestake are over the experimental ones for all values of sin 2 (20), but as this 
parameter tends to 1 the fit of the values turns better and the rates of Homestake can be well 
explained within 3.5<r for sin 2 (20) > 0.75, considering the rates with or without recoil-electron spec- 
trum[|. This is also verified in Super-Kamiokande for sin 2 (20) > 0.5. For the gallium experiments all 
rates are checked for any sin 2 (20) with 3.5cr, but small values of this parameter give a better agree- 
ment, so we can check GALLEX and SAGE within an experimental 2a error for sin 2 (20) < 0.25. 

1 Note that Homestake, GALLEX and SAGE are not affected individualy by the recoil spectrum as Kamiokande 
and Super-Kamiokande are, but we are dealing with the data from the plots that fits the five experiments together. 
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Table 4: Suppression rates using the averaged solar neutrino survival probabilities, with 99% C.L., 
considering maximal suppression of 7 Be neutrinos for Homestake and gallium experiments and min- 
imal suppression of pp neutrinos for gallium. 



sin 2 (20) 


Homestake 


Gallium 


rates only 





0.44 ±0.16 


0.67 ±0.14 


0.25 


0.38 ±0.11 


0.60 ±0.10 


0.50 


0.36 ±0.08 


0.51 ±0.06 


0.75 


0.35 ±0.04 


0.44 ±0.02 


1 


0.34 ±0.02 


0.42 ±0.01 


rates + recoil-electron spectrum 





0.43 ±0.17 


0.66 ±0.14 


0.25 


0.38 ±0.12 


0.59 ±0.11 


0.50 


0.36 ±0.09 


0.49 ±0.07 


0.75 


0.34 ±0.05 


0.42 ±0.02 


1 


0.33 ±0.02 


0.42 ±0.01 



For Kamiokande the data are explained for all values of sin 2 (20) within 3.5<r, but for sin 2 (20) = 1 it 
is checked with 2a. 

A better fit for Homestake is obtained considering the maximum possible suppression of 7 Be 
neutrino^ (this suppression is limited by eq.|l5]); in special for sin 2 (20) > 0.5 this agreement is 
found whitin 3a error. Unfortunatelly this fact leads to a decrease in the GALLEX and SAGE 
rates. To explain these rates with the maximum suppression of 7 Be neutrinos, we need to assume a 
(Pu e u e )p P > 0.60. In this case the fits are obtained, within 3a error, for any values of sin 2 (20). See 
table 4 for more details. 

So a better agreement for the five experiments is found, with 99% C.L., for sin 2 (20) > 0.50. 
The allowed regions that fit the observed rates of the five experiments with 99% C.L. are shown in 
the figure 5. In this figure are plotted the calculated rates using the maximum suppression of 7 Be 
neutrinos and (P Ue u e )pp — 0.60. From the figure we see that the allowed regions of Sm^i x sin 2 (2o;) 
grows when sin 2 (20) rises. For sin 2 (20) : 
10- n eV 2 ; for sin 2 (20) 



0, the range of the parameters are, 0.61 < sin (2a;) < 1 



and 8m\ x > 5.0 



0.50, 0.48 



< sin 2 (2o;) < 



1 and Sml-L > 4.8 



10- LL eV 



and for sin 2 (20) = 1, 0.29 < sin 2 (2o;) 



< 1 and 8m 2 2l > 3.1 



10 eV . The x analysis for all 



plots in figure 5 are shown in table 5 and we can observe again that the inclusion of the recoil-electron 
spectrum gives allways smaller values of x 2 - 



2 The total suppression of the 7 Be neutrino only occurs for sin 2 (20) = (see figure 2). 
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Table 5: Best fit points for the Sm^ x sin 2 (2u;) plots shown on figure 5. Column 1 is for the seasonal 
effect with the rates only and column 2 is with the inclusion of the recoil-electron spectrum. 



sin 2 (20) 


X 1 

A.min 


sin 2 (2cu) 


5m 2 2l (eV' 2 ) 


Column 1 





8.89 


0.62 


8.32 . 10- 11 


0.25 


9.64 


0.96 


7.60 . 10- 10 


0.50 


10.5 


0.86 


7.48 . 10- 10 


0.75 


11.5 


0.91 


3.47 . 10- 10 


1 


13.9 


0.72 


3.47 . 10- 10 


Column 2 





4.86 


0.63 


8.90 . 10- 11 


0.25 


5.40 


0.56 


8.90 . 10- 11 


0.50 


6.13 


0.48 


8.90 . 10- 11 


0.75 


6.98 


0.42 


2.75 . 10- 10 


1 


7.56 


0.44 


4.42 . 10- 10 



5 Conclusions 

The recent Super-Kamiokande data on solar neutrinos brought new constraints on the previously 
allowed solutions and have restricted them. In this work we show that the vacuum oscillations still 
give a good explanation to the solar neutrino problem and that three neutrino solutions gives better 
results. In three neutrino scenario the allowed regions for Sm^ and sin 2 (2a;) parameters that fits all 
experimental data increase as sin 2 (20) grows. 

The calculated suppression rates including the recoil-electron spectrum leads to a reduction of 
approximately 30% in the mass regions, but it does not affect the fitting of the data because the 
averaged probabilities are almost the same as those obtained considering the rates only. An analysis 
with and without the recoil-electron spectrum shows that the inclusion of this spectrum checks better 
the experimental data. 

A better fit for the experiments are obtained, with 99% C.L., if we consider the maximum 
suppression of 7 Be neutrinos and a suppression of pp neutrinos from to 40% as well a suppression of 
about 60% of 8 B neutrinos. In this case the best agreement with the experimental data is obtained 
for sin 2 (20) > 0.50. The inclusion of the recoil-electron spectrum gives a better x 2 fit than the 
calculated with rates only, for any value of sin 2 (20). 
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99% C.L. in vacuum, for different values of sin 2 (20). The full lines represent the best fit for P VeVe 
calculated with the rates only (column 1 on figure 3) and the dashed lines the same fits with rates 
and recoil-electron spectrum (column 2 on figure 3). 

Figure 5 - Mass regions versus mixing angles for the BP98 SSM with 99% C.L. in vacuum, with 
seasonal effect, for different values of sin 2 (20) considering maximal suppression of 7 Be neutrinos. 
The first column is for the rates only and the second column is for the rates with recoil-electron 
spectrum. The black dot in each plot represents the best fit. 
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